Abstract. An automated search for star clusters close to 
the Galactic plane (\b\ < 5°) was carried out on the Point 
Source Catalogue of the DENIS survey. 44% of the Galac- 
tic plane have been observed and calibrated. The method 
allowed to retrieve 22 known star clusters and to iden- 
tify two new ones, not published yet although previously 
presented in the 2MASS web site as embedded clusters in 
HII regions. Extinction in the field and in front of the clus- 
ters are estimated using a model of population synthesis. 
We present the method and give the properties of these 
clusters. 
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1. Introduction 

The highly concentrated dust in the plane of the Milky 
Way makes difficult the observation in the visible range of 
embedded objects and most distant objects hidden behind 
the Galactic plane. However, the extinction decreases at 
longer wavelengths: in the K-band, it is about ten times 
smaller than in the visible. Thus the recent Near-Infrared 
surveys provide suitable data for systematic search for new 
objects projected on the Galactic plane. Us ing the Two 
Micron A ll Sky Survey (hereafter 2MASS, jBkrutskie ei] 
al., 1997] ) i n the J (1.25 fim), H (1.65 fim) and K s (2.17 
/j,m) bands, Dutra & Bica (2000) carried out a systematic 
search for new star clusters in a 5° x 5° field centered 
close to the Galactic plane and listed 58 IR star clus- 
ters or candidates. They also investigated selected regions 
with evidence of star formation and found 42 star clusters 
or candidates ( putra fc Bica, 2001[ ). Hurt ct al. (2000D 



scrcndipitously discovered two globu lar clusters near the 
Galactic plane. Vaughn ct al. (2001 ) undertook a search 
for galaxies by eye on J and K s -band images from the 
Deep-Ne ar In frared Survey of t h e Sou thern Sky (hereafter 



1999) and found 13 star 



DENIS, |Epchtein et al., 1997 , 
clusters (Rousseau, private communication). 

Using the data in I (0.82 /im), J, and mostly K s bands 
from DENIS, we performed an automated search for star 
clusters in a ±5° band around the Galactic plane. In 
Sect. |^ we describe our method and detection criteria. 
The basic properties of the new star clusters are given in 
Sect. H 

2. Method 

DENIS observations are strips 12' wide in right ascension, 
and 30° spanned in declination. A strip is made of 180 
individual 12' x 12' images, with 2' overlaps between con- 
secutive images. The images are processed at the Paris 
Data Analysis Center, and point source catalogues are ex- 
tracted. Over the 3317 strips calibrated at that time, we 
selected the strips which contain sources with Galactic lat- 
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itude \b\ < 5°, that is 1219 strips. They all together cover 
44% of the entire Galactic plane, with Galactic longitudes 
ranging between —160° and 40° (82% of the strips are 
calibrated within this range). 

The two criteria used to search for star clusters are 
based on the density of stars and the integrated flux in 
the K s band. This band allows to probe more extinguished 
regions than the / and J bands. Each image of the strip 
is divided in 2' x 2' frames. The number of stars N and 
total flux F in the K s band are computed in each frame, 
using the corresponding calibrated catalogue. The mean 
N and F, and the dispersion around the mean o~n and 
up of these two parameters are computed for the entire 
image. A 2' x 2' window is then displaced with 0.4' steps 
in declination and right ascension on the image. If the 
density of stars N and the total flux F in the window are 
4 xctjv and 4 xap above the mean N and F, the program 
indicates the eventuality of a detection. The threshold has 
been defined with the faintest globular cluster discovered 
with 2MASS, GC01. 

However, the star cluster GC01 is not easily distin- 
guishable only from the catalogue, and these two criteria 
only are not sufficient to perform an automatic detection. 
Most of the stars in the cluster have J-K s > 3.5, which is 
not surprising for a globular cluster dominated by red gi- 
ants in a high extinction region. With the condition J-K s 
> 3.5, the cluster can now be detected in an automatic 
way. A difficulty remains in the determination of the J-K s 
limit due to the unknown color of an eventual cluster and 
to the variable extinction from frame to frame. We chose to 
consider the middle of the color range of the stars, J-K mi( j, 
in the 12' x 12' image as a tentative limit. The efficiency 
of the detection being very dependent on the selection of 
the limit J-K m id considered, we performed the 4 xer/v and 
4 xerp test 3 times, over stars redder than J-K m jd - 0.5, 
J-K,„id, and J-K m id + 0.5. In the case of GC01, the de- 
termined limit is J-K m id = 3. The detection program fails 
with the condition J-K m ;d > 2.5 or 3, but is successful for 
J-K m id > 3.5. The variation of the J-K s limit increases the 
efficiency of detection, but also allows to detect younger 
star clusters than this globular cluster. 
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Wilson (1991) and [Battinclli (199l| ) described auto- 
matic procedures for identification of star associations. 
They are based on agglomerative algorithms. Whereas 
Wilson (1991) determined a search radius, Battinclli 
(199l|)J algorithm, called Path-Linkage-Criterion tech 
niquc, assigns two stars to the same group if it is possible 
to go from one to the other jumping from star to star in 
steps smaller than a determined scale-length. It allows to 
detect compact and filamentary groups. Our method is not 
optimized to identify associations of young stars. It uses a 
color criterion in order to enhance the contrast with back- 
ground stars, combined together with a density criterion 
and a total luminosity criterion. Applied to Near-Infrared 
data, our technique allows to retrieve either old globular 
clusters or young star clusters embedded in HII regions. 

3. Results 

Already known star clusters in the studied region detected 
with this method are listed in Table 1. Most of the non- 
detected clusters are faint and do not appear clearly in 
the DENIS images. NGC 6544 and NGC 6553 were not 
detected whereas they are bright globular clusters. The 
first one is effectively not revealed in the catalogue prob- 
ably because the bright part in the center has been ex- 
tracted as a single bright source. The second one is in 
between two strips and only partially visible in one strip, 
whereas the adjacent strip is not available. Two star clus- 
ters not yet published have been identified. However, these 
clusters have been previously presented in the "Picture of 
the week" 2MASS web site, as embedded clusters in HII 
regions Gum 25 and W40. Properties of these clusters, 
hereafter Gum 25 cluster and W40 cluster, are given in 
Table 2. Images in the K s -band are shown in Fig. |]. (K S ,J- 
K s ) color-magnitude diagrams are plotted in Figs. |p and 
||a for Gum 25 cluster and W40 cluster respectively. Field 
stars are plotted with dots, cluster stars with pluses. 

3.1. Determination of cluster characteristics 

While the DENIS data alone do not allow to determine the 
definitive membership of a star to a cluster (kinematics or 
distances would be needed), one may estimate a probable 
membership from the position on the sky and estimate the 
contamination by field stars. From the density around the 
clusters, the expected contamination in both clusters is 5 
and 6 field stars respectively. The cluster size is estimated 
visually from the image and catalogue, as well as the stel- 
lar content. We have attempted to estimate the reddening 
of the clusters and in the field by fitting an extinction 
variable along the line of sight. The Besangon model of 



population synthesis (Haywood et al., 1997) is used to re 



produce the stellar content in the fields. The extinction 
on the pathsight is fitted in order to reproduce the stel- 
lar distribution in the (K S ,J-K S ) plane. It is done using a 
distribution of diffuse extinction which is produced by a 



Table 1. Known clusters in the studied region detected 
with the method described in Sect. E| *C: star clusters. 
OC: open clusters. GC: globular clusters. 



name 




a 


(J2000) 


5 




1(°) 


b(°) 




LDN1654 a 


06 


59 


42 


-07 


46 


29 


220.79 


-1 


71 


*C 


Pismis 2 


08 


17 


55 


-41 


40 


18 


258.86 


—3 


34 


*C 


NGC 2660 


08 


42 


18 


-47 


09 


00 


265.85 


—3 


03 


OC 


Pismis 12 


09 


19 


54 


-45 


08 


00 


268.64 


+3 


19 


*c 


Westerlund 2 


10 


24 


02 


-57 


45 


30 


284.27 


-0 


33 


OC 


NGC 4337 


12 


23 


54 


-58 


08 


00 


299.29 


+4 


54 


*c 


NGC 5927 


15 


28 


44 


-50 


40 


22 


326.70 


+4 


79 


GC 


NGC 5946 


15 


35 


29 


—50 


39 


35 


327.58 


+4 


19 


GC 


Lynga 7 


16 


11 


03 


-55 


18 


52 


328.77 


-2 


79 


GC 


NGC 6256 


16 


59 


33 


-37 


07 


17 


347.79 


+3 


30 


GC 


Tonantzintla 2 


17 


36 


11 


-38 


33 


13 


350.80 


-3 


42 


GC 


Terzan 5 


17 


48 


05 


-24 


46 


48 


3.84 


+1 


69 


GC 


NGC 6440 


17 


48 


53 


-20 


21 


34 


7.73 


+3 


80 


GC 


NGC 6441 


17 


50 


13 


-37 


03 


04 


353.53 


-5 


01 


GC 


2MASS GC01 


18 


08 


22 


-19 


49 


47 


10.47 


+0 


10 


GC 


2MASS GC02 


18 


09 


36 


-20 


26 


44 


10.07 


-0 


45 


GC 


Terzan 12 b 


18 


12 


16 


-22 


44 


31 


8.36 


-2 


10 


GC 


NGC 6712 


18 


53 


04 


-08 


42 


21 


25.35 


-4 


30 


GC 


NGC 6749 


19 


05 


15 


+01 


54 


03 


36.20 


-2 


20 


GC 


NGC 6760 


19 


11 


12 


+01 


01 


50 


36.10 


-3 


92 


GC 


NGC 6530 


18 


04 


48 


-24 


20 


00 


6.14 


-1 


38 


*C 



a (Hodapp (1994J ) 

b also designed as ES0522SC1 

Table 2. Gum 25 cluster and W40 cluster characteristics. 
Both are embedded in indicated HII regions. 



Gum 25 cluster 



W40 cluster 



a,5 (J2000) 

l,b 

size 

A v 

HII region 

— distance 

- Vlsr d 



09 02 11,-48 49 14 

269.27,-1.53 

1.4' x 1.2" 

9 mag 

Gum 25 a 

1.8 kpc 

7 km s _1 



18 31 25,-02 05 02 

28.80,+3.50 

1.9' x 1.3' 

17 mag 

W40 b 

600 pc 

0.6 km s" 1 



a also designed as RCW40 

b also desi gned as Sh2-64, RCW174,[L 89bl28. 790+0.346 

c Gum 25: |Avedisova fe Palous (1989| ) W40: |Wu et al. (1996| ) 

d Gum 25: |Brand et al. (1987|), W40T[Lockma,n (198S|) 



double-exponential disc of interstellar material, plus dis- 
crete clouds as needed on the line of sight. The overall 
extinction at position (Z, b) and distance d writes: 

f d 

A v (l,b,d) = / A dif x x e^^dr + V] Cl(l, b, d) 

where hi and hh are the extinction scale length and scale 
height, Cl(l,b,d) the extinction due to the intervening 
cloud in (I, b) direction at distance d, (R, z) the galactic 
coordinates. The scale length and scale height are taken 
as 4.5 kpc and 140 pc, following the expected distribution 
of the dust. Ada is the differential diffuse extinction at the 
solar position in magnitude per kpc. This value is often 
taken between 0.5 and 1. This range of value works well 
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Fig. 1. DENIS K s -band images of the new star clusters. Left: Gum 25 cluster. Right: W40 cluster. The size of the 
images is 5.8' x 5.8'. Top = North, left = East. 
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Fig. 2. Color-magnitude diagram (K S ,J-K S ) for Gum 25 cluster, a: Observed diagram. Field stars are represented with 
dots, cluster stars with pluses. Photometric errors are 0.05 at K s = 12 and 0.15 at K s — 13. The dashed line shows 
the J completeness limit, b: Simulated diagram (see text), c: Observed diagram for cluster stars only, superimposed 
with Padova isochrones of different ages. 



(a) 



+ 
+ 



. :.- i- >«;/•• • t • 
+ ■ . a%\ ■: 



2 3 
J-K 



(b) 



2 3 
J-K 



(c) 



+ 
+ 



4Myr nOOMyr/ 1Gyr 




2 3 
J-K 



Fig. 3. Same as Fig. || for W40 cluster. Photometric errors are 0.2 at K s — 11 and 0.3 at K s = 13. 



at high and intermediate latitudes but often failed in the 
plane. It so need to be adjusted on each line of sight. 

The distribution of extinction described above gives 
a univoque relation between distance and extinction in a 
given field. However observations sometimes show rapid 
variations of extinction from one line of sight to another 



especially in dense clouds. This can be accounted for by as- 
suming a dispersion of extinction for a ny given star (each 
star being on a different line of sight). Lada et al. (1994 ) 
found in the cloud IC5146 variations of extinction which 
can be modeled with dispersion proportional to the ex- 
tinction: a = 0.73 + 0.40 x Ay in a field of 1.5' on a side. 
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Thoraval et al. (1997) find a lower value: a = 0.25 x Ay. 
In adjusting our model of extinction we have tentatively 
used either a constant extinction at a given distance or 
variations proportional to the extinction with the above 
formula. 

3.2. Gum 25 cluster 

In this field it has been easy to describe the extinction 
with a simple model of diffuse extinction distributed over 
the line of sight with no high density cloud. The (K S ,J-K S ) 
diagram of field stars is well reproduced with a differential 
diffuse extinction at the solar position Am = 2.5 mag 
kpc -1 (Fig. I^d). With this model, the extinction reaches 
about 10 magnitudes at 10 kpc. 



Using theoretical Padova isochrones (Bertelli et al 



1994; Girardi et al., 1996), we try to estimate the red- 



dening of the cluster. However it appears difficult to fit a 
well defined isochrone on the color-magnitude diagram of 
this cluster due to the large dispersion in colour. The age 
cannot be constrained either. This dispersion cannot be 
interpreted by field star contamination alone and must be 
due to a strongly varying extinction inside or in front of 
the cluster. We estimate the mean extinction to be Ay=9 
mag. Fig. ||c shows the cluster (K S ,J-K S ) diagram together 
with the isochrones at a distance of 1.8 kpc, which is the 
distance of the Gum 25 HII region, and with a visual ex- 
tinction of 9 magnitude. This is higher than the value 
proposed by 2MASS on their web page, A v —3.6 mag, a 
value which we propose as a lower limit for less extincted 
stars in the cluster. 

3.3. W40 duster 

The distribution of extinction in this field has been more 
difficult to estimate. The simple extinction model from 
the interstellar matter disc does not allow to explain why 
there arc so few main sequence stars in this field. A high 
extinction cloud at 600 pc from us, which is the distance 
of the W40 ionized region, allows however to reproduce 
the number of main sequence stars if this cloud has a vi- 
sual extinction of 11 magnitudes. In this case the other 
background stars are reddened to J-K > 2. The (K S ,J-K S ) 
diagram of simulated stars is given in Fig. ||b. Main se- 
quence stars are at distances less than 2 kpc, while giants 
appear at 4-8 kpc. At a distance of 6 kpc in the direction 
of W40 cluster, the stars are above the diffuse extinction 
disc. Thus the diffuse extinction is not sufficient to explain 
the large dispersion among field giants. In order to repro- 
duce the distribution in J-K s of the giants as well as the 
observed dispersion, we adopt an internal dispersion on 
the extinction at a given distance of a — 0.25 x Ay. The 
distribution of extinction along the line of sight is defined 
by Aiif = 3 mag kpc -1 , and 11 mag extinction due to a 
cloud at 600 pc. With this model, the extinction reaches 
about 21 magnitudes at 10 kpc. 



In the cluster, the isochrone fitting of the main se- 
quence is obtained with the cluster being at the same 
distance as the main cloud d = 600 pc and an extinc- 
tion in the cluster Ay = 17 mag, larger than the cloud 
which reddens the field. It can be interpreted as the cen- 
tral dense part of a larger cloud. Fig. ||c shows the main 
sequence fitting of the cluster with Padova isochrones of 
4 Myr, 100 Myr and 1 Gyr, visual extinction of 17 mag 
and distance of 600 pc. It shows that the cluster is rather 
young. 

4. Conclusion 

We elaborated a program to perform an automated search 
of star clusters from the catalogues of extracted sources 
from the DENIS survey. We concentrated on the band 
around the Galactic plane (\b\ < 5°), where 44% of the 
data have already been calibrated. Most of the already 
known star clusters visible on the DENIS images in the 
probed region have been detected. Two star clusters, not 
published yet, have been identified. Both clusters are em- 
bedded in their associated HII regions. Extinction in front 
of the clusters is estimated. Uncertainties remain large on 
this determination mainly due to questionnable member- 
ship and to the photometric errors at faint magnitudes. 
Deeper and more accurate photometry and proper mo- 
tions would allow to assert the cluster characteristics. 

Although we have missed known star clusters and 
probably not yet detected clusters, our method applied 
to Near-Infrared data allows to search for old star clus- 
ters and embedded star clusters in a systematic way, with 
a lower efficiency than when looking at the image, but 
much more rapidly. We plan to run again this program 
once the complete region around the Galactic plane will 
be calibrated. 
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